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The small- a; behavior of structure functions in the satu- 
ration region is determined by the non-linear generalization 
of the BFKL equation. I suggest the effective field theory 
for the small- a; evolution which solves formally this equation. 
The result is the 2-1-1 functional integral for the structure 
functions at small x. 

PACS numbers: 12.38. Bx, ll.lO.Jj, 11.55.Jy 

The great success of pQCD in describing the be- 
havior of structure functions of deep inelastic scattering 
(DIS) can be traced back to the fact that the depen- 
dence is governed by DGLAP evolution equations which 
have two remarkable properties: they are linear equa- 
tions, and the evolution at high is purely perturbative 
(the non-perturbative physics enters the game only when 
we lower the normalization point /i^ down to the typical 
hadronic scale ~ IGeV). The higher-order terms of per- 
turbative expansion for both the coefficient functions and 
the anomalous dimensions of the light-cone operators lie 
in the same framework of linear evolution and lead to the 
corrections ^ cts,ct1 etc. 

The situation for the small-x DIS is more complicated. 
In the leading logarithmic approximation (LLA) the 
small-x asymptotics is described by the BFKL pomeron 
101 . It is possible to reformulate the BFKL equation as 
an evolution equation where the relevant operators are 
Wilson lines - infinite gauge links |^ The evolution 
of the two- Wilson-line operator ( "color dipole" ) with re- 
spect to the slope of Wilson lines reproduces the BFKL 
equation. 

Unfortunately, the theoretical status of the BFKL evo- 
lution is not as clear as the DGLAP one (for the review, 
see Refs. Q). The biggest problem is the lack of uni- 
tarity: the power behavior of the BFKL cross section 
violates the Froissart bound and therefore, in order get 
the true asymptotics at small x, we must go beyond the 
LLA. At this step, we face a new problem. In the DGLAP 
case, the sub-leading logaritms follow the same general 
pattern of linear DGLAP equation and the problem is 



^ At high energies the particles move so fast that their trajec- 
tories can be approximated by straight lines coUinear to their 
velocities. The proper degrees of freedom for the fast particles 
moving along the straight lines are the (infinite) gauge factors 
ordered along the straight line M. 



purely technical: calculating the loop corrections to the 
kernels. In the case of small- a; evolution there are also 
as corrections to the BFKL kernel but, on the top of 
that, there are the unitarity corrections which lie outside 
the framework of the BFKL equation. At small Us and x, 
these corrections seem to dominate over the NLO BFKL 
ones Q. 

Another problem with the BFKL evolution is infrared 
instability. We can safely apply pQCD to the small-x DIS 
if the characteristic transverse momenta of the gluons 
fcx in the gluon ladder are large. For the the first few 
diagrams, one can check by explicit calculation that the 
characteristic fc^ are ~ Q^. However, as x decreases, it 
turns out that the characteristic transverse momenta in 
the middle of the gluon ladder drift to Aqcd making the 
application of pQCD questionable. 

Recently, an idea has emerged that these two difficul- 
ties may cancel each other out. Consider the DIS from 
the heavy nuclei where the large density sets the satura- 
tion scale Qs which effectively cuts the integration 
over fcj^ even at relatively low energy. As we shall see be- 
low, the small-x evolution in this case is non-linear which 
leads to the growth of the saturation scale with energy, 
see the discussion in Refs. [0-^. It is natural to assume 
that even for the DIS from the nucleon where there is no 
saturation at low energies, the saturation scale at suffi- 
ciently small X may be generated by the non-linear evo- 
lution itself. Indeed, the parton recombination described 
by the non-linear evolution must balance at some point 
the effects of parton splitting so the partons will reach 
the state of the saturation. In this high-density regime 
the coupling constant is small but the characteristic fields 
are large, making a perfect case for the application of the 
semiclassical QCD methods mi^,|l^. The high-density 
regime of QCD can serve as a bridge between the domain 
of pQCD and the "real" non-perturbative QCD regime 
governed by the physics of confinement. 

In this paper I suggest the effective field theory which 
describes the small-x evolution in the saturation region. 
First, let me remind the OPE for high-energy amplitudes 
derived in Consider the amplitude of forward 7*7*- 

scattering at small xb ~ In the target frame, the 
virtual photon splits into qq pair which approaches the 
nucleon at high speed. Due to the high speed the classi- 
cal trajectories of the quarks are straight lines collinear 
to the momentum of the incoming photon q. The corre- 
sponding operator expansion switched between nucleon 
states has the form [|j: 
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ld^xe^'^-^{p\T{j^ix)MO)}\p) 

d^x^I^,{x^){p\TT{U{x^)U\0)}\p), 



(1) 



where Ifj,u{xi_) is a certain numerical function of the 
transverse separation of quarks Xi_ and virtuahty of 
the photon = —q^. The relevant operators C/(f/^) 
are gauge factors ordered along the classical trajectories 
which are almost light-like lines stretching from minus to 
plus infinity: 



\ J — OO 
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where e is collinear to q and z± is the transverse position 
of the Wilson line. 

It turns out that the small- a; behavior of structure 
functions is governed by the evolution of these opera- 
tors with respect to the deviation of the Wilson lines 
from the light cone; this deviation serves as a kind of 
"renormalization point" for these operators. At infinite 
energy, the vector e is light-like and the corresponding 
matrix elements of the operators (|^) have a logarithmic 
divergence in longitudinal momenta. To regularize it , 
we consider operators corresponding to large but finite 
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velocity and take e(; — ei + C,e2 where ei — {q— 
and e2 = p are the lightlike vectors close to the directions 
of the colliding particles. Now, instead of studying the 
energy-dependence of the physical amplitude we must in- 
vestigate the dependence of the operators (Q) on the slope 
Large energies mean small C, and we can sum up log- 
arithms of C instead of logarithms of s (At present, we 
can do it only in the leading logarithmic approximation 
(LLA) as <C 1, In ^ ~ 1). The equation governing 
the dependense of J7 on C has the form ||^ |^ 



C,—U{x^,y^) 
dC, 



(3) 
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dz I 



+ U{z±,y±_) -U{x^,y±_) +U{x^, z^)U{z_i_,y±_)^ 

where U{x±,y±) ee ^ Ti{U{xj_)W {y±)} - 1. The first 
three linear terms in braces in the r.h.s. of eq. (j^)) 
reproduces the BFKL pomeron ||] while the quadratic 



■^The first non-linear equation for parton densities is known 
since 1983 as the GLR equation (it was conjectured in Ref. [Q] 
and proved in the double- log limit in Ref. [^). The full LLA 
X result was first derived in Ref. by t he above method. 
After that, it was reobtained in Ref. in the framework 
of the dipole model , in Ref. by direct summation 

of relevant Feynman diagrams, and in Refs. by the 

semiclassical methods. 






FIG. 1. BFKL evolution in terms of Wilson-line operators 
(denoted by dotted lines). 



term will give us the three-pomeron vertex |T8|. The 
solution of the linearized evolution equation is especially 
simple in the case of zero momentum transfer (e.g. for 
the total cross section of small-x DIS): 



{p\U^-^-{x^,0)\p) 



dv 
2^ 
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C] and m? is ei- 



where uj{u) = 2Nc^[-Rexl:{\ 
L/ 

-). The sketch of linear evolution is 



ther or (in LLA, we cannot distinquish between 



Q^s In Tw E^nd Us In ■ 
presented in Fig. 1. The starting point of the evolution 
is the slope collinear to the momentum of the incoming 
photon q (C = xb) and it is convenient to stop the evo- 
lution at a certain intermediate point C,o = ^ where 
So ^ ^ In <C 1 . The first of these conditions 

means that sq is still high from the viewpoint of low- 
energy nucleon physics while the second condition means 
that So is sufficiently small from the viewpoint of high- 
energy physics (so one can neglect the BFKL logs). The 
matrix element of the double- Wilson-line operator at this 
slope is a phenomenological input for the BFKL evolu- 
tion (just as the structure function at low serves as the 
input for ordinary DGLAP evolution). At large s the in- 
tegral over v is dominated by the vicinity oi v = Q which 

gives the familiar BFKL asymptotics (T*°' ~ x^"^ " 

Unlike the linear evolution, the general picture is very 
complicated since the number of operators U and in- 
creases after each evolution. At the time being, it is not 
known how to solve the non-linear evolution equation in 
an explicit form. It is possible, however, to write down 
the solution of the non-linear equation (||) in the form of 
a functional integral over the double set of the variables, 
<rj=i,2(^±,??) = t°'^f=i.2i^A^,'n) belonging to the Lie alge- 
bra of the SU(3) color group and ^i=i,2{z^,'n) belonging 
to the group itself: 

C/'W(a;_L)«)i7t''^(y±) = (5) 
/ Z?^i(z, v)Dc,2{z, ii)Dniiz, ii)Dn2{z, 77) 

"'f2l,2(')0) = l 

n\{x^,7^A)UTn2{x^,iiA)<)^nl{v±.VA)ul^''ni{v±.iiA) 

expJ r drjjd^zl- ^ <,t{z.v)d'Ul\{z,i^)h,{z,rfX 
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where tt = and {n^iiy = 2Tr{rrjtf2}. Going to 

the the variables tt = we see that Eq. (||) is a phase- 
space functional integral for the non-local Hamiltonian 

H{7:uTr2,ni,n2) = (6) 

dx^dy^7:nx±)ix^\^[dl{n\n2y'']^\y±)7rl{y^) 
P± P± 

where |x| is an eigenstate of the coordinate operator nor- 
niaUzed according to = S'^^^x — y), see e.g Ref. 

pl| . The rapidity 77 serves as a Euclidean "time" for this 
evolution. 

We shall demonstrate that the perturbative expansion 
of the functional integral (m reproduces the evolution of 
the color dipole U{x±) ® U'iy^) in the LLA. To get the 
perturbative series, we substitute fl{x±, rf) — e~'^3'l>{x±,n}-^ 



Tri,2{VA)=0 



ni=i,2-DQ(z, r])D(f)i{z, r]) 



01,2(')o)=O 



X y 

X e-^'^'^i'^^'''-^) expJ dri I (f-z\- ^ ^"(2,77) 
\ dri J 

- ^q"(z,r;)<^2^2,77)92(^e*9'^l(^.'?)[/;^Og-^902(.,r,)^"''j| 

Next, we can represent the r.h.s. of Eq. (|^) in the form 

1,2(»?a)=0 

n,^i,2DQ{z,r])D(j),{z,7]) (8) 

<Ai,2(r;o)=0 

[nA,Vo]xUS''{VO, Va}x {VA, Vq}vUI'^" [VO,VA]y 

X expJ f drjfd'^z - i ^ <:,°-{z,r])d^tj)i{z,ri) 



- rj)d'{[riA, mhU^"{Vo, VAj^Y^^i^, v) 

where we introduced the notations 

r 1 rri «S f" ^ 1 . , is f' " 02 (a: i , J?) 

(9) 

Let us now expand the r.h.s. of Eq. in powers of g. 
The first nontrivial term in this expansion is 



(10) 



/■7ri,2(i;A)=0 

as i U.,=i^2DQ{z,ri)D(j>i{z,ri) 

J 0l,2(r;o)=O 

(^My±,VA}ul'^° - ul^'>My±,VA)] 

(^Mx±,Vo)UT(t>2{x±,VA)) (EUl"^^' 

UT ® {My±,VA)ul'^<'My±,VA)) 

dr, I d'z,^{z,7^)4{z,T,)^(u^o{z^)y' 
r 

dr^ d^z J2 'it{z,l)d^k{^,il) 



X exp 



The propagators for this phase-space functional integral 
are 

(0^(x,,7)^}(2/,r/)) = »<5,^.J-&((a;|l|y))% - 77'), 
(0^(x,r;)05(y,V)) =0, (^f(x,r;K}(y,r;'))=0 (11) 
With these propagators, the r.h.s of Eq. ([l^) reduces 

to 



-as(j]A - fyo) 



X (xl^d^U'^" ^\yY'' ~ eUS°t''^Ul'^°{x\^{d^U''«) 
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(12) 



which coincides with the Eq. (B17) from Ref. |^]. Taking 
trace over the color dipole indices one reproduces the Eq. 

. Similarly, it can be demonstrated that further terms 
of the expansion of eq. (H) in powers of g repoduce the 
subsequent iterations of the non-linear equation (^. 

The intergral over tt variables can be easily performed 
resulting in: 

U'^^{xi:)®U^'^^{yi_) (13) 
rf2i.2(r;A)=0 

Dni{z,rj)Dn2iz,Ti) nl{x±,TjA) 

I Ol,2(lo) = l 

X t/"" {x^)n2{x^,'nA) ® r!^(y±, t^a);/^"" (y±)r!i(2/±, 7?a) 
X exp I 5- / df] I d^ 



Vo 



d''(n\{z,ri)Ut^'>ni{z,ri) 



ah 



X 92 (^zr!l(z,r;)jii(z,,7))"92(^zr!t(z, 77)572(2,77))'} 

Note that the action of this effective field theory is local. 
This functional integral for the small- a; evolution of the 
Wilson-line operators is the main result of the paper. 



3 




FIG. 2. Propagation of the color dipole through the nu- 
cleus. 



In the case of large nuclei it is possible to write initial 
conditions for the small-x evolution using the McLerran- 
Venugolalan model. The nuclear matrix element of the 
two- Wilson-line operator ( "color dipole" ) is given by the 
Glauber formula, see Fig. 2. 



cPzj_{A\TrU{x_L + z^)U{zi_)\A) 



(14) 



Here Li, = — 6^ is the propagation length of the 

dipole (located at the impact parameter h) through the 
nucleus, p 



4/37ri?3 



4(iV2 



is the nuclear density, and 



(15) 



The Eq. (^ is derived under the assumption that the 
characteristic size of the dipole (the "saturation scale " ) 
is smaller than the size of the nucleon ^ In this case, 
the quarks propagating along the straight light-like lines 
^ interact by the instantaneous (in the light-cone time 
x_|_) potential 



d'pA 



(2^)2 2p4 



(16) 



where toq 



2^ is the IR cutoff. Q It is worth noting 
that the factor —1 in the parenthesis in the l.h.s. comes 



This assumption is certainly true at yl oo. For real 
nuclei, one should find the saturation scale Qa from the final 
result for the matrix element of the color dipole between the 
nuclear states, and verify that Qs S> IGeV 

^As we mentioned above, the energy soq should be high 
enough so we can replace the slope pi + C,op2 by pi in the 
non-logarithmical expressions. 



from the diagrams with the two gluons attached to the 
same nucleon and the same Wilson line. Taking into 
account the color factors, one obtains the Eq. m£) with 
xbG(xb, p^ = ^I^) = ^lii2;J^/mo, see Ref. |2^] . 

Similarly to Eq. ([l3|), it is possible to represent this 
result as a functional integral over a variable A(a;j^,^) G 
5f7(3): 



(17) 

DA{y, I) A{x + z, Lb) A^z, Lb) 



d'z{A\UT+M'"\A) 

A'(Li„y)=0 
A(0,a) = l 

X exp[^l^'^'dlJd'y{A{l,y)A'{l,y)y 

X {-d'+ml)^{A{l,y)A'{l,y)y} 

where A' = ^A. Extra J/''" (a;) (U^^'^ix)) lead to extra 
A{x,Lb) [A^x,Lb)) in the pre-exponent. 

The final formula for the martix element of the color 
dipole operator at small xs is obtained by combining the 
functional integrals (O) and (p^): 



d^z{A\U''^''{x + z) ® C/t''^ {z)\A) 

pA'{Lb.,v±)=0 r(li,-2{riA,v)=0 

dH / DAil,y)Dn,iy,rj) 

JA{Q,y^) = l Jnt,2{rio.y) = l 

X DVl2{y,-q) Vl\{x + z,r]A)A{Lb,x + z) 

(g) VL2{x + z,riA)^l{z,riA)A\Lb, z)VLi(z,riA) 



X cxp -j^ - 
1 f"-^ 



1 



dl / d^y{iA{l,y)A'{l,y)Y 



'^g^P Jo 

X (ml~^)^tA{l,y^)A'{l,y)y 
dvjd^yd^(in\iy,ri)n,{y,r^)y 

P(nliy,ri)A^Lb,y)niiy,v)] 

X ^{inliy,r^)n2iy,v)] 




(18) 



The gluon structure function in the LLA is propor- 
tional to the matrix element of the dipole operator 
xbG{xb,P^ = xl^) = -^(A|Trf/;'-^(xi)C/r(0)|A)^o 
the numerical calculation of the functional integral (flq ) 
should give the nuclear structure functions at small x. 
This would be complementary to the approximate solu- 
tions of Refs. ||ll|,|l^,^2|jl9|j2^@ since it could give the 
structure functions not only in the asymptotic black- 
body limit, but also in the intermediate region defining 
the saturation scale Qs- 

It should be mentioned that our formula ( p^ ) gives 
the evolution of the color dipole only in the LLA. In the 
case of large nucleus we have an additional parameter 
yl 3> 1 so our LLA approximation based on the non-linear 
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equation has a window a'^A^^^ 1, a^lnxs 1 
where it is justified even at moderately small xb- In the 
case of nucleon, our as{Qs) 1, Q;s((5s) Inxs ~ 1 ap- 
proximation should be justified a posteriori after check- 
ing that the saturation does occur at sufficiently small 
Xb- If the saturation takes place at such low x that 
as{Qs) lii-XB S> 1, our LLA breaks down and we need to 
take into account the non-fan diagrams such as t-channel 
loops formed by BFKL pomerons. However, the non- 
linear equation (||) leads to the result for the structure 
function which does not violate unitarity (see the discus- 
sion in Refs. ]ri|- p^ , p7|p5| ) and therefore we should not 
expect the large discrepancy between the unitary LLA 
result and the exact amplitude at present energies. 

Acknowledgements. The author is grateful to Y.V. 
Kovchegov, E.M. Levin and L. McLerran for valuable 
discussions. This work was supported by the US Depart- 
ment of Energy under contract DE-AC05-84ER40150. 



[17: 

[is; 

[19 

[20 

[21 

[22; 

[23' 
[24; 

[25' 
[26; 



184 (1994); Z. Phys. C64, 631 (1994). 
M.A. Braun, Eur. Phys. J. ClQ(2Qm) 337. 
J. Bartels and M. Wusthoff, Z. Phys. C66, 157 (1995). 
E. lancu and L. McLerran, Preprint SACLAY-TOl-026 



( |hep-ph/0103032| ); 

H. Weigert , Preprint 



ph/0004044) 



NORDITA-2000-34-HE ( pp 



I. Balitsky and V.M. Braun, Nucl. Phys. B311, 541 
(1989). 

B.Z Kopeliovich, I.L. Lapidus, and ALB. Zamolodchikov, 
JETP Letters 33(1981) 612. 



A.M. Mueller, Nucl. 
Yu.V. Kovchegov 
B529(1998) 451. 
A.M. Mueller, Nucl. 
E. Levin and K . 
(|hcp-ph/0012167| ); 



ph/0101275) 



Phys. B335 (1990) 115; 

and A.M. Mueller, Nucl. Phys. 

Phys. B437 (1995) 107. 
Tuchin, Preprint TAUP-2659-2000 
Preprint TAUP-2664-2001 (|hep 



[1] V.S. Fadin, E.A. Kuraev, and L.N. Lipatov, Phys. Lett. B 
60(1975), 50; I. Balitsky and L.N. Lipatov, Sou. Journ. 
Nucl. Phys. 28(1978), 822 . 
L Balitsky, Nucl. Phys.B 463 (1996) 99. 
L.N. Lipatov, Phys. Reports 286 (1997) 131; 
I. Balitsky, "High-Energy QCD and Wilson Lmes" , in 
the Boris loffe Festschrift "Handbook on QCD" , ed. M. 



Shifman (W orld Scientific, Singapore, 2001), e-print [lep- 
ph/0101042||. 



O. Nachtmann, Annals Phys. 209, 436 (1991); J.C. 
Collins and R.K. Ellis, Nucl. Phys.B 360 (1991) 3. 
V.S. Fadin and L.N. Lipatov, Phys. Lett. B 429(1998) 
127; G. Carnici and M. Ciafaloni, Phys. Lett. B 
430(1998) 349 . 

Y.V. Kovchegov, A.M. Mueller, Phys. Lett B 439(1998) 
428 . 

L.V. Gribov, E.M. Levin, and M.G. Ryskin, Phys. Kept. 
100(1983) 1 . 

A.H. Mueller and J.W. Qiu Nucl. Phys. B268(1986)427. 
L. McLerran and R. Venugopalan, Phys. Rev. D49 (1994) 
2233; Phys. Rev. D49 (1994) 3352; Phys. Rev. D50 
(1994) 2225. 

A.H. Mueller, Nucl. Phys. B558 (1999) 285. 

Yu.V. Kovchegov, Phys. Rev. D6G, 034008 (1999); Phys. 

Rev. D61,074018 (2000). 

E. Levin and K. Tuchin, Nucl. Phys. BB573, 833 (2000). 
J. Jalilian-Marian, A. Kovner, L. McLerran, and H. 
Weigert, Nucl. Phys.B 504 (1997) 415; Phys. Rev. 
D59(1999) 014014; Phys. Rev. D59(1999) 034007; 



L Balitsky Phys. Rev. Lett. 81(1998) 2024, [jhep 

ph/9808215t , Phys. Rev. D60, 014020 (1999). 

A.H. Mueller, Nucl. Phys. B415(1994) 373: A.H. Mueller 

and Bimal Patel, Nucl. Phys. B425(1994) 471. 

N.N. Nikolaev and B.C. Zakharov, Phys. Lett. B 332, 



5 



